Bosonization for disordered and chaotic systems 
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Using a supersymmetry formalism, we reduce exactly the problem of electron motion in an external 
potential to a new supermatrix model valid at all distances. All approximate nonlinear a models 
obtained previously for disordered systems can be derived from our exact model using a coarse- 
graining procedure. As an example, we consider a model for a smooth disorder and demonstrate 
that using our approach does not lead to a "mode-locking" problem. As a new application, we 
consider scattering on strong impurities for which the Born approximation cannot be used. Our 
method provides a new calculational scheme for disordered and chaotic systems. 
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The supersymmetry method has proven to be a use- 
ful method for solving different problems for disordered 
systems 0. Within this technique, for the problem of 
weakly disordered conductors, one derives a supermatrix 
nonlinear a model which is valid for distances exceed- 
ing the transport mean free path i tr and is applicable for 
a description of the diffusive motion of electrons. The 
main object in this approach is an 8x8 supermatrix field 
Q(r) that depends on the coordinate r. The free energy 
functional for space variations of Q(r) has the form 



f[q] = y J Str [ D{ yQ) 2 + 2 *( w + iS ) A Q] dr > 



(1) 



where D is the classical diffusion coefficient, v is the den- 
sity of states at the Fermi surface, and w is the frequency. 
The matrix A = diag(l, — 1) distinguishes between ma- 
trix elements corresponding to advanced and retarded 
Green functions and Str means supertrace. The super- 
matrix Q(r) satisfies the relation Q 2 (r) = 1. 

Recently, a lot of attention has been paid to the study 
of ballistic electron motion and relatedproblems of quan- 
tum chaos (for a review see, e.g., Ref.|2(). Equation £[J is 
not applicable in this case and one needs a generalization 
of the nonlinear a model. It is expected that the proper 
model can be applied to clean systems which are chaotic 
in the classical limit. Treating disordered and chaotic 
systems in the same manner is of great advantage since 
they share common properties such as the universality of 
the level statistics, although the origin of randomness is 
totally different. 

In these problems, the traditional derivation of the 
nonlinear a model based on the disorder averaging is of 
no use and one needs more general derivation. Several 
attempts to generalize the a model and include shorter 
distances have been undertaken in the past several years 
El IEL El • Although this study has lead to consider- 
able progress in understanding the ballistic motion, the 
goal of the complete description has not been achieved. 
We mention that the most general quasiclassical form of 
the ballistic u model has been written for a smooth disor- 
der potential in a recent work Q . It is valid for distances 



down to the wavelength Ap. However, it seems that even 
shorter distances are also important because in the ab- 
sence of a disorder, the approach of Ref. (as well as all 
others) does not allow one to derive a so-called "regular- 
izer" of Refs.JEEl that is believed to originate from the 
quantum diffraction and is absolutely necessary to obtain 
the Ehrenfest time. 

In this Letter, we suggest a scheme of derivation of 
the a model that is comple tely different from all those of 
Refs. 0, 0, 0, IE IE S IE IS 113 In contrast to the previ- 
ous derivations, we are able to reduce the initial electron 
problem to a model containing an integral over superma- 
trices Q( r iP) exactly. This model is a generalization of 
the o~ model to all distances and is written for any given 
potential. It can be simplified by integrating out shorter 
distances and this is how semiclassical ballistic a models 
and the diffusive model can be derived. 

The method we use is some kind of the bosonization 
that maps the original fermionic system onto a bosonic 
one, which has been used for a wide variety of physical 
systems O ur approach is also known as the geo- 

metrical quantization or the coherent-state method and 
is very close to that in Ref.0] (a similar trick has been 
used in a proof of the universality for random matrices 
H3)- At the same time, the supersymmetry allows us to 
do considerably more and obtain a comparatively simple 
explicit form for the bosonized action. 

Using 8-component supervectors ip(r) we write the 
generating function Z[J] for the electron systems as Q 
Z[J] = J exp(-£)X>(</>, i>) with 



C = -i I rb(r)Hj (r,r')V>(r')drdr', 

lo + iS 



(2) 

A J -J(r,r'), (3) 



Hj = 5(r-r')[H-E 



where H is the Hamiltonian of the electron system and 
J(r,r') is a source field. The variables e and u> stand 
for the electron energy and the frequency of the external 
field, respectively. We write the Hamiltonian as H — 
p 2 /2m—eF+U(r), where ep is the Fermi energy and U(r) 
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is a potential that may include a regular part as well as 
an impurity potential. An extension to more complicated 
systems including magnetic or spin-orbit interactions is 
straightforward. 

Now we want to express the original generating func- 
tion Z[J] in a "bosonized" form. We start our derivation 
with the identity 

1= J e -ij Str [(l/2)Q(ry)-i ) (v)^(v')]A(v' ,r)drdr' v q va ^ 

where Q and A are arbitrary supermatrices with the same 
symmetry as the product ?/>(r)?/>(r'). This identity is ev- 

I 



All functions in the above equations are written using 
the Wigner representation A(R, p) = J dpA(r, r')e _lpp , 
where R = (r + r')/2 and p = r — r', and we use the 
notation x = (r, p) for a point in the phase space. The 
products of two supermatrices A(x) * B(x) are defined 
by the Moyal formula 

A(x) * B(x) = A(x)eW*X*'V»-™»)B(ic). (7) 

Calculation using the Wigner representation and the star 
product * is known as the Weyl symbol calculus 0] and 
is convenient for quasiclassical expansions. 

The generating function Z[J] is expressed in terms of 
the integral over the variables Q and A. Remarkably, 
due to the supersymmetry one can calculate the integral 
over A exactly. This can be done writing the superma- 
trix Q(x) as <5(x) = T(x) * q(x) * T(x), where q(x) is a 
diagonal matrix and T(x)*T(x) = 1. Changing the vari- 
ables A(x) — > T(x) * A(x) * T(x) we express the integral 
B as a function of the variable q(x) only. We can proceed 
further making one more change of the variables of the 
integration A(x) — > q~ 1/2 (x) * A(x) * q~ 1/2 (x). The Ja- 
cobian arising when changing the variables is unity due 
to the supersymmetry. Then we find the new variable 
A does not couple to Q. Because of the supersymmetry 
the final integral over A is exactly equal to unity and we 
obtain the generating function Z[J\ = j cxp(— F[Q])VQ, 
where 

F[Q] = l - J Str [Hj(x) * Q(x) - i In Q(x)} dx. (8) 

This is the main result of the present work. We see that 
the initial electron problem which is written in terms of 
the supervector ip( r ) is reduced exactly to a functional 
integral over the supermatrices Q(x) depending on the 



ident because the integral over the variable A gives the 
delta function S(Q/2 — ipijj). The factor 1/2 is written 
to keep along the notations of Ref. Q . Substituting this 
expression into Z[J] we integrate over the supervectors 
ip(r) and shift the variable A as A — > A sp + A, where 
A sp = Hj — iQ^ 1 ■ Then we obtain 

Z[J] = J exp{-| j Str[Hj(x)*Q(x)]dx\B[Q]VQ, 

(5) 



(6) 



coordinates x in the phase space. Since ^(r) are su- 
pervectors and not fermionic vectors, we may call the 
procedure leading to Eq.JHJ "superbosonization" . 

We emphasize that the supersymmetry is used in this 
formulation for writing the generating function in terms 
of an integral over matrices but not for averaging. No av- 
eraging over impurities has been performed in the deriva- 
tion of the functional F[Q) and the derivation itself is 
just a sequence of identical transformations. This con- 
trasts the traditional supersymmetry method pj where 
the supersymmetry was used for averaging over disorder 
in the beginning of the calculations. The variables Q(x) 
in F are different from the variables of integration used 
in other related works and the representation (jSJ, being 
as exact as the original representation, is simply more 
convenient for making approximations. Its validity is by 
no means related to a form of the disorder distribution. 

Now we want to show how the formalism works in dif- 
ferent situations. First, let us consider a system with a 
smooth potential such that one may expect a quasiclas- 
sical behavior. We demonstrate how one can come to the 
model of Ref.0 starting from Eq.JSJ. We do not discuss 
here clean billiards and how one should average over the 
spectrum. This is an important and complicated problem 
but it is beyond the scope of the present Letter. Instead, 
we consider, as in Refs.0, 0j an infinite system with a 
smooth disorder. 

The functional F[Q] is quite complicated and we find 
first its minimum in Q(x) at J = 0. The minimum is 
achieved at qo(x) = iHj= Q = g(x), where g(x) is the 
one particle Green function of the Hamiltonian without 
sources (multiplied by i). If the potential U{r) in the 
Hamiltonian is smooth, the dependence of g(r, r') on (r + 
r') /2 is weaker than on r — r' and we can write with good 



B[Q]= / exp 



Str{ln[iQ- x (x) - A(x)\ - M(x) * Q(x)}dx)VA 
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accuracy 



7T J n J 



£(p) + «A 



A. 



(9) 



where £(p) = p 2 /2m — £p. Writing this equation we 
assumed that the spectrum is continuous, which is def- 
initely true for an infinite metallic sample. In a finite 
system, the use of g(x) in the form of Eq.(J5J may be 
justified after a certain averaging over impurities. 

In order to compute the contribution of fluctuations 
near the saddle point we write Q(x) in Eq.© as Q(x) = 
T(x) * q(x) * r(x), where q(x) is a block-diagonal super- 
matrix (commuting with A) and T * T = 1. If J = 0, 
fluctuations of the supermatrix g(x) near qo do not give 
any contribution and its average is equal to the exact one 
particle Green function g(x). This is because fluctuations 
of T do not contribute to the one particle Green functions 
1] and one can simply put T = 1 calculating this quan- 
tity. Then we immediately come to the conclusion that 
(<?(x)) = <?(x), which means that the integration over the 
"massive" modes (fluctuations of q) does not give any 
correction. 

In contrast, fluctuations of T(x) are important and 
the main contribution comes from matrices slowly vary- 
ing in the phase space. We perform first the integration 
over supermatrices q by making a cumulant expansion in 
gradients of T. The main contribution comes from the 
first cumulant. Putting q(x) = <?(x) in Eq.(8j, we can 
easily integrate over |p| because only g(x) should be in- 
tegrated, whereas the momentum p in the supermatrices 
T(x) should be taken at the Fermi surface, p = pfh, 
n 2 = 1. Expanding the star products in gradients, which 
corresponds to the quasiclassical approximation, we re- 
duce F to the form 

F = T / Stl ' { Af " (r) ^ FnVr ~ 1 V r C/(r)V n ] T n (r) 



+i 



lu + iS 



A - J„(r) 



Qn(r) >drdn 



(10) 



where <2 n (r) = T n (r)AT n (r). This functional agrees with 
the corresponding functional derived from an equation 
for quasiclassical Green functions 0]. Equation (|10|l is 
applicable for all distances exceeding the wavelength Af- 
One can find a detailed investigation of properties of this 
functional and its reduction to Eq.Q in Ref.0- 

It is clear from the derivation that the momenta p are 
pinned in the supermatrices T(x) to the Fermi surface 
due to the singular form of the supermatrix g(x) = <?(x) 
near the Fermi surface. In the traditional approach based 
on the Hubbard-Stratonovich or color-flavor transforma- 
tion, the matrix q(x) is smooth at the Fermi surface and 
nothing suppresses fluctuations perpendicular to the en- 
ergy shell (the "mode-locking" problem). We see that 
the mode-locking problem is related merely to an in- 
convenient choice of the variables of the integration in 



the traditional approach and this problem does not exist 
here (it does not exist in the diagrammatic approach 
and the approaches based on the quasiclassical equations 

S0). 

In the absence of the random potential U(r) one has 
only terms linear in gradients and no regularizer such as 
that suggested in Refs. 0,01 El appears in Ea. lfTUI) . which 
would mean that classical orbits do not mix. However, 
we are sure that a more careful consideration of quantum 
effects leading to the mixing of the classical motion is 
possible within the model (JSJ. 

Now let us consider a disordered metal with a small 
concentration of strong short-ranged impurities, such 
that the self-consistent Born approximation(SCBA) can- 
not be used. This problem cannot be considered by the 
traditional a model scheme of Ref. [jj because the saddle 
point approximation used in this approach is equivalent 
to the SCBA. Moreover, we are not aware of any diagram- 
matic calculation of the weak localization correction for 
this model, although the effect of str ong impurities was 
considered in different situations 0, 

Our derivation of the reduced a model consists of three 
steps: the coarse-graining procedure, the quasiclassical 
approximation, and the impurity averaging. First, we 
integrate out short distances of the order of Af and re- 
duce the initial model to an effective one describing fluc- 
tuations at longer distances. After this coarse-graining 
procedure, we can use the quasiclassical approximation. 
At the same time, the impurities are self-averaged over 
the integrated length scale. It is important to notice that 
this method of the averaging is different from the tradi- 
tional method where the impurity is averaged by consid- 
ering the ensemble of systems. Since all of the variables 
are slowly varying in space, we can simply average the 
terms of the coarse-grained functional by changing the 
observation point, which is equivalent to the averaging 
over the positions of the impurities. 

In order to perform the coarse-graining procedure, we 
rewrite the supermatrix Q as Q(x) = V(x) * Q(x) * V"(x), 
where V = V^ 1 describes slow modes and Q stands for 
fast modes. Substituting this representation into Eq.lJBJ 
we write Q(x) = g(x) + SQ(x) and expand F in SQ up to 
the second order. Taking into account fast modes [fluctu- 
ations of <3(x)] is necessary to obtain the correct form of 
the collision term. This is done by performing Gaussian 
integrals over SQ. Alternatively, we could use the con- 
traction rules for Q which can be obtained by comparing 
Z[J] written in terms of ip and Z[J) of Q. This allows us 
to take into account higher order effects in a systematic 
way. As a result, we obtain 

F = T / StrA( 3 dx+ ^ / Str #*Qdx 

~/str (g*V*[H*V]) 2 dx, (11) 
where the Q matrix has been reduced to Q(x) = V(x) * 
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g(x) * V"(x). We neglected the source term for simplicity. 

Since the variables V(x) are slow we can use the qua- 
siclassical approximation. This means that we introduce 
the quasiclassical Green function <? n (r) = 7r _1 J g(x)<i£ as 
Eq.@ and expand the star product. The Green function 
g(x) still includes the impurity potential and we intro- 
duce the function X(x) as g(x) = go(p) * [1 + iX(x) * 
5o(p)]i where go(p) is the one particle Green function 
in the absence of impurities. Writing the impurity po- 
tential as a sum over the positions r of the impurities 
U(r) — J2a u ( r ~ r i)' wnere u ( r ) is the potential of a 
single impurity, we average over r a . A standard consid- 
eration (see, e.g., [l^) leads for a small impurity con- 
centration n; to the following result in the momentum 
representation 

M ^-i-£mI»«>-*- (12) 

where / is the exact amplitude of scattering on a sin- 
gle impurity. Then we find for the Green function 
(g(p + q/2, p - q/2)) imp = *d(q)/(£(p)+iA/2T), and for 
the quasiclassical Green function, g n (r) = A. We intro- 
duced the scattering time r using the optical theorem 
Im/(n,n) = %vA J dn'|/(n, n')| 2 = {2n l r)- 1 A. Other 
contributions such as the average of a product of X can 
be considered in a similar way and are expressed in terms 
of the scattering amplitude /. 

We can write the coarse-grained functional using 
the scattering amplitude /(p, p') and the free particle 
Green function go{p) instead of using the impurity poten- 
tial U(r) and the full Green function g(x). Introducing 
the quasiclassical Green function, the momentum p is re- 
stricted to the Fermi surface. We write the supermatrix 
Via the form V n (r) = F(r)[l-iP n (r)-(l/2)p2(r) + - • •], 
PA + AP = and expand Eq.^TTJ in P. We obtain 

F = T / Stl '{ A ^( r >F nV r^n(r) + yAQ n (r) 

+w(n,n') [P n (r) - P n ,(r)] 2 JcZrdndn', (13) 

where iu(n,n') = 2Trv \f(n, n')| 2 is the probability of 
scattering from the state with n to the state with n', and 
Qn(r) = V n (r)AV n (r). In Eq.JTSJ), we keep in the colli- 
sion (last) term second powers of P only. When averaging 
over the impurities, it is difficult to keep all powers of P 
and this is the reason why we do not write higher order 
terms. Equation (|13fl is not sufficient for the complete de- 
scription of the ballistic regime where the collision term 
is, generally speaking, a complicated function of P. How- 
ever, as long as we are interested in the diffusive regime, 
Eq. (jT3|l is enough and we can come to Eq. Q by a Gaus- 
sian integration over the supermatrices P. The diffusion 
coefficient D in Eq.£Q) is given by D = VpT tr /d (d is the 



dimensionality), where the transport time r tr is defined as 
t^ 1 =n; / iu(n,n')(l-nn')dn', andQ(r) = V(r)AV(r). 

In conclusion, we developed a new approach to disor- 
dered and chaotic systems. The most general superma- 
trix model obtained in Eq.© is exact and valid at all 
distances. It can be simplified by using a coarse-graining 
procedure for both smooth disorder and strong impuri- 
ties, which leads to ballistic and diffusive nonlinear a 
models. We believe that our method is a proper route to 
a more complete description of the quantum chaos. 
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